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ABSTRACT 

There has been recent speculation (Davies & King 2005) that the cores of intermediate-mass 
stars stripped of their envelopes by tidal interaction with the supermassive black hole in the 
Galactic centre could form a population observationally similar to the so-called Sgr A* cluster 
or 'S' stars, which have close eccentric orbits around the hole. We model the evolution of such 
stars, and show that the more luminous end of the population may indeed appear similar to 
young B stars within the observational limits of the Galactic Centre region. Whether some or 
all of these cluster stars can be accounted for in this manner depends strongly on the assumed 
IMF of the loss cone stars and the scattering rate. If most of the observed stars are in fact scat- 
tered from the Galactic Centre inner cusp region itself then the population of ~ 20 to current 
observational limits may be reproduced. However, this only works if the local relaxation time 
is small and relies on the cusp stars themselves being young, i.e. it is dependent on some star 
formation being possible in the central few parsecs. Conversely, we obtain a possible con- 
straint on the tidal stripping rate of 'normal' -IMF stars if there are not to be red stars visible 
in the Sgr A* cluster. 
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1 INTRODUCTION 

In recent years, our ability to make observations of the region 
around the central black hole of our Galaxy has significantly im- 
proved. What we see there, however, remains somewhat of a para- 
dox. It is becoming apparent that the stars we observe closest to the 
centre are predominantly hot, blue stars - even though relatively 
unobscured observations can only be made in the H and K bands 
which are strongly biased towards picking up those stars which are 
redder - and that their appearance suggests that they are young and 
massive. Genzel et al. (2003) observe two disks of hot stars within 
the central 1" - 10" (0.04 - 0.4 pc), which appear to contain ab- 
normally high numbers of Wolf-Rayet (WR) and luminous blue 
variable (LBV) stars (Paumard et al. 2001), both products of the 
later stages of massive stellar evolution. There is some evidence 
that both this region and the Galactic Centre area in general have 
initial mass functions (IMFs) which are either biased towards the 
high-mass end or else have an abnormally high lower mass limit 
with respect to a standard Salpeter-like IMF (Nayakshin & Sun- 
yaev 2005, Stolte et al. 2005). However, the extreme tidal forces 
in the Galactic centre region make the normal formation of stars in 
situ a difficult task. Either these apparently massive stars must have 
formed in a non-standard manner or else they must have migrated 
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inwards in a timescale short enough to accommodate the lifespan 
of an O star (Hansen & Milosavljec 2003), a few 10 7 years at most. 

The Sgr A* cluster is located even closer to the central black 
hole than the hot star populations described above - it contains 
the stars with the smallest known distances of approach to Sgr A* 
(Ghez et al. 2003). Unlike the disk stars, their orbits are isotrop- 
ically distributed and highly eccentric (e ~ 0.3 - 0.9), with peri- 
bothra as low as 0.5 mpc. These stars are also blue. K-band spec- 
troscopy (Ghez et al. 2003; Eisenhauer et al. 2005) suggests that 
they are young O and/or B stars. The formation problems associ- 
ated with the massive stellar disks are even worse for the Sgr A* 
cluster. If these stars formed in situ, their formation was subject to 
even stronger tidal forces. Similarly, any inward migration is likely 
to have been even more strongly impeded. It has been suggested 
that they are the remains of a shredded open cluster drawn inwards 
by dynamical friction (Kim & Morris 2003), perhaps containing an 
intermediate mass black hole to keep it tightly bound enough to get 
in close enough to deliver them to their current location (McMillan 
& Portegies Zwart 2003). However, these solutions would lead to 
spatial distributions of stars different from those observed (Ghez 
et al. 2005). Perhaps the most promising solution is that the high- 
velocity stars (HVSs) observed to be heading away from the Galac- 
tic Centre with velocities of hundreds of kms~' (e.g. Brown et al 
2005) share a common origin with the Sgr A* cluster, via the split- 
ting of binary systems which encounter the black hole. In this see- 
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nario, one component of the binary is ejected to become a HVS 
and the other is captured (Gould & Quillen 2003; Perets, Hopman 
& Alexander 2006). 

One possible solution to the problem is that the Sgr A* cluster 
stars are not young stars at all - that they merely look like them. 
As proposed by Davies & King (2005), another way to make blue 
stars in the Galactic centre region is to strip away the outer layers 
of giants, exposing their hot cores. On many of the orbits on which 
the Sgr A* cluster stars currently reside, a RGB or AGB star would 
overfill its Roche lobe by several times at the point of closest ap- 
proach to the black hole. Given that another of the well-known mys- 
teries regarding the central region is the lack of red giants therein 
(Genzel et al. 1996; but see also Jimenez et al. 2006), this suggests 
that any giant star which ends up in the region after scattering onto 
an orbit passing close to Sgr A* quickly loses some or all of its en- 
velope through tidal interactions with the black hole. The exposed 
core of such a star will appear blue, and it is possible under the 
restricted observing conditions imposed by the Galactic centre en- 
vironment that it may resemble an O or B star for the portion of its 
stripped lifetime that it is visible. 

This formation route for the Sgr A* cluster stars also accounts 
for their observed tightly-bound, eccentric orbits. Although the 
scattered star initially approaches the black hole on a near-parabolic 
orbit, the energy required to unbind the envelope is of a similar or- 
der to the orbital energy loss required to produce orbits similar to 
those observed. This energy requirement also puts constraints on 
which progenitor stars can feasibly reproduce the observed Sgr A* 
cluster population. Input stars must be able to survive the strip- 
ping process without being completely shredded. Those for which 
the stripping of the entire envelope does not provide enough en- 
ergy loss to produce a closely-bound orbit will end up on nearly- 
parabolic orbits in which they are only observable close to pericen- 
tre - a negligible fraction of their lifetimes. Input stars must also 
have large radii to end up on orbits similar to those observed. The 
initial distance of closest approach to the black hole corresponds to 
the periastron distance in the post-stripping orbit. At periastron dis- 
tances similar to those observed, a main-sequence star is not large 
enough to fill its instantaneous tidal lobe. Such stars have to pass 
far closer to the central black hole to undergo tidal mass loss. This, 
combined with their less centrally-condensed structures, suggests 
that main-sequence stars are probably shredded completely once 
they have filled their tidal lobes. Envelope stripping in the central 
region, albeit by stellar collisions, has also been invoked by Alexan- 
der (1999) to explain the distribution of stars in the inner 0.05 pc. 

There are further observational constraints that these stripped 
giants must obey if they are to be a plausible explanation for the 
Sgr A* cluster. The cluster stars all appear to be on or near the 
main sequence; therefore we require that the observable population 
of stripped stars must also primarily reside in this area of the HR 
diagram. Since K-band visibility increases rapidly with decreasing 
T e ff, this is (if one ignores the potential effect of low-luminosity 
AGN irradiation on their atmospheres; Jimenez et al. 2006) a strong 
constraint on how red the stars may be. However, a population bluer 
than the main sequence is simply invisible in the K-band to current 
observational limits and is therefore only constrained by the re- 
quirement that the scattering event rate (Wang & Merritt 2004) and 
the mass of unobservable remnants (Goodman & Paczynski 2005) 
be realistic. The stripped stars must also display a similar spec- 
trum in the K-band to those observed. This requires, for instance, 
a composition including a non-negligible surface hydrogen abun- 
dance, and an absence of Wolf-Rayet-like emission lines. Lastly, 
they must also be capable of reasonably rapid rotation; for exam- 



ple, the star S2 (S0-2) has a vsim value of around 220km s . Of 
these, the event rate is the most pressing problem. In this paper we 
consider whether those constraints can be met by running simula- 
tions of rapidly-stripped stars from different input populations and 
comparing them with the observations. 



2 SIMULATIONS 

The process of tidal stripping is a complex one, and it is far from 
certain that the apparent balance between envelope binding energy 
and orbital energy translates into the possibility of a real and ef- 
ficient transfer of energy between the two. To study this, we take 
a two-pronged approach. Firstly, the tidal mass loss itself and its 
effect on the orbit is being studied by means of detailed SPH sim- 
ulations (Dale et al., in preparation). Secondly, the evolution of 
stripped stars with simple assumptions about the stripping process 
over a range of input parameters is calculated in order that their 
properties may be compared with those of the Sgr A* cluster stars. 
This is the subject of the current paper. In doing so we have as- 
sumed the stripping process itself and the uncertainties relating to 
it (whether, for example, there is really an efficient conversion be- 
tween orbital and envelope binding energy, and whether the stellar 
core is able to survive the process) to be a 'black box' to be ex- 
plained or disproved by the SPH simulations. 

2.1 Code 

We use the Eggleton stellar evolution code (Eggleton 1971, Pols 
et al. 1998), including the convective overshooting prescription of 
Pols et al. (1998) with 8 0V = 0.12, which is roughly equivalent to 
an overshooting distance of 0.2H p , where H p is the pressure scale 
height. This leads in general to relatively large model cores. All of 
the runs in this paper were carried out at a metallicity of Z = 0.02 
as appropriate for the Galactic Centre region (Najarro et al. 2004; 
Carr, Sellgren & Balachandran 2000). A slightly higher or lower 
metallicity is unlikely to make a strong difference. Wind mass loss 
is included according to the prescription of de Jager et al. (1988), 
but is generally small for most of the masses we have considered. 
Whilst we can run the code following binary evolution including 
mass loss to a companion (in this case the black hole), we here 
treat our star as a single star, and a constant high rate of mass loss 
to it at the appropriate age to simulate the stripping event. 

The rapid nature of the mass loss from the stars we are con- 
sidering presents a problem for simulation. Essentially, the entire 
envelope may be removed in a single orbit of the black hole; this 
equates to a mass-loss rate which may be over a solar mass a year, 
which is difficult to apply in any realistic way to a stellar evolution 
code. What we can do, however, is apply mass loss at the fastest 
rate which the code will sustain. This is an approximate approach, 
as equilibrium is assumed throughout, but gives a general idea of 
how these stars will behave. The mass-loss rate used varies slightly 
from star to star - since we use in each case the largest rate the 
code will support without breaking down - but is generally between 
10~' and 1O~ 3 M yr _1 . In all cases the entire mass loss event lasts 
significantly less long than would be a typical timestep for normal 
evolution of that star at that evolutionary stage. From the evolution 
of other stars which are or have been affected by rapid mass loss 
- the primaries in interacting binary systems, subdwarf O and B 
stars and massive Wolf-Rayet stars - it is possible to infer that the 
most likely course for our stripped stars is a rapid journey to the 
blue part of the HR diagram followed by a relatively stable period 
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Figure 1. Typical evolutionary tracks for intermediate mass stripped stars 
pre- and post-stripping. Shown in each case are a 3.25 Mq (lower track) and 
an 8.85M .; model (upper track) Pre- stripping tracks are shown in black, 
the rapid mass loss and subsequent adjustment phases are indicated by the 
grey lines and post-stripping tracks are shown in blue (8.85 Mq) and red 
(3.25 M r ). The stage at which stripping occurs and the amount of envelope 
remaining are given on each panel. 

of helium-burning on the He main sequence, as can be seen in Fig. 
1. 

2.2 Model Grid 

A number of variables will determine the evolution of a star after 
its envelope has been removed. 

2.2.1 Initial Mass 

There is no absolute upper or lower mass limit on stars which in- 
teract with Sgr A*. However, the range of plausible pre-stripping 
masses for at least the visible cluster stars is limited at the lower 
end by stars which are unlikely to become luminous enough to be 
observable in the Galactic centre region, and at the upper end by 
stars which are likely to have a Wolf Rayet-like appearance when 
stripped and are thus relatively easily distinguishable from B stars. 
We limit our detailed simulation input mass range to between 1 and 
50Mq. In this range we run model grids at 21 log-spaced masses, 
and interpolate as required between the resulting tracks. However, 
we consider stars over the whole range of potential IMFs by us- 
ing the synthetic stellar evolution formulae of Tout et al. (1997) for 
lower masses. For these stars it is highly unlikely that we can create 



a blue core which is as luminous as those observed - a rapid tran- 
sition to white dwarf status is the most likely outcome of stripping. 
Furthermore, most stars below lM e have sufficiently long evolu- 
tionary timescales that at the current age of the Universe they will 
still be on the main sequence, and are therefore ruled out as progen- 
itors for at least the currently-observed Sgr A* cluster stars. Whilst 
it is possible that the Sgr A* cluster stars are stripped stars in the 
upper region of or even above the mass range we simulate, the short 
lifetimes of very massive stars (a few Myr) are a strong impediment 
to getting them in to the Galactic centre region in time. It is still pos- 
sible, however, that they could have been scattered from the disks 
of stars at ~ 0.4pc, which are massive and may have been formed 
in situ (Nayakshin & Sunyaev 2005). Normally intermediate-mass 
stars are strongly favoured by the IMF over more massive stars, 
but there is evidence that in the Galactic Centre region the IMF is 
significantly top-heavy or may have a high lower limit (Stolte et 
al. 2005). However in this case the original (albeit not the current) 
masses of the Sgr A* cluster stars would be greater than that usu- 
ally inferred in the case that they are B stars! In this scenario the 
'Paradox of Youth' (Ghez et al. 2003) is not resolved; at the least it 
is a requirement that the massive star disks formed in their current 
location. Such an origin may also require that the relaxation time in 
the central parsec is shorter than the life of a ~ 6M0 star, i.e. ~ 10 s 
years. 

The simulations are ended either when the star has become 
a white dwarf, when degenerate carbon burning ignites or when 
central carbon burning ends for the more massive models. 



2.2.2 Evolutionary Stage at Stripping 

Since tidal stripping can occur rapidly when a giant is scattered 
onto a orbit which passes close to the central black hole, the evo- 
lutionary point at which the stripping occurs does not necessarily 
correspond to the stages in the star's life when it is beginning to 
expand, as it would for Roche lobe overflow (RLOF) in a binary 
system. Instead, this violent RLOF-like event may occur at any 
point in the star's life when its radius is above the instantaneous 
tidal radius on the scattered orbit, 

where p is the periastron and M* and Mbh are the masses of 
the scattered star and the black hole respectively. However, anal- 
ogously to RLOF, there is also the possibility that stars which are 
scattered onto near-parabolic orbits when their radii are too small 
for the stripping process to take place (as is likely for main se- 
quence stars, which should form the bulk of the scattered popula- 
tion) may be stripped on a subsequent orbit of the black hole when 
they have evolved into larger-radius objects. This requires that the 
orbital time be less than the evolutionary timescale of the star - in 
the case of our best-matching models a period of a few 10 6 years is 
probably the maximum for these orbits, corresponding to an initial 
tidal interaction on first pass close to the black hole which extracts 
some 10 4f) ergs of energy from the orbit. Depending on the relax- 
ation time in the region near Sgr A* and the resulting period, some 
of these stars are also likely to be scattered back off these long or- 
bits before they can be stripped. 

To allow for a distribution of scattering times, we simulate for 
each of our stellar models the onset of rapid mass loss at the RGB 
base, middle and tip, on the horizontal branch and at the base and 
middle of the AGB. Stars which approach close enough to tidally 
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overflow whilst they are still on the main sequence are likely to be 
completely disrupted. 

2.2.3 Amount of Envelope Remaining 

Without detailed calculations of the effect the tidal forces have on 
the star (currently in progress) it is hard to say exactly how much 
of the envelope will be stripped. For the vast majority of possible 
orbits the answer will be 'nearly all' but the remaining amount can 
strongly affect how blue the star becomes, if it becomes blue at all, 
and also its remaining lifetime after the stripping event. In order to 
get an idea of how much of the envelope needs to be removed be- 
fore a blue star is obtained, we run models at each of the above grid 
points which are left with 0, 1, 5, 10, 20 and 50 percent of their en- 
velopes. The latter few situations are unlikely - requiring scattering 
onto a grazing orbit and therefore relatively mild mass loss. Typi- 
cal evolutionary tracks are shown in Fig. 1. Generally stars with 
more than about 5 percent remaining envelope spend a significant 
amount of their subsequent evolution as red giants. 

The amount of envelope left strongly affects the orbit the 
stripped star ends up on, as it limits the amount of energy which 
can be extracted by envelope removal. It is possible that different 
Sgr A* cluster stars have different amounts of remaining envelope. 
It is also notable that since the smaller the mass of the stripped core, 
the smaller its orbital energy (in an equivalent orbit), the criterion 
that the binding energy of the removed envelope be comparable 
to the orbital energy a typical Sgr A* cluster star would have if it 
were at the remaining core mass is fulfilled by a wide range of ini- 
tial masses, provided the star is large enough to overflow its tidal 
lobe in that orbit. This is illustrated in Fig. 2 for stellar models at 
various masses. 



3 INPUT POPULATION 

Simulations which calculate the disruption rate of stars by black 
holes frequently restrict the stellar input distribution to a single 
mass and radius (e.g. Wang & Merritt 2004). However, that ap- 
proach is clearly of little use here, since we wish to simulate a pop- 
ulation which is likely to arise from a relatively small segment (the 
high-mass, large-radius end) of a diverse input population. Input 
masses will initially be distributed according to an IMF which, at 
least for stars originating from outside the immediate Galactic Cen- 
tre region, is likely to be Salpeter-like (e.g. Kroupa 2001). However, 
the probability of a star ending up on an orbit which passes close to 
Sgr A* also depends on its mass and radius, as does its final orbit. 
If the star is not stripped on its first pass of the black hole, whether 
it can be stripped on a subsequent pass depends on its future radius 
evolution. 

3.1 Scattering and event rate 

The tidal disruption loss cone of a galaxy is the set of orbits for 
which stars pass sufficiently close to the central black hole for tidal 
disruption to take place. This can only happen once per star so, as 
well as having to be populated in the first place, the loss cone must 
be continually refilled if the number of tidal disruptions is not to 
drop to zero with time (e.g. Merritt & Wang 2005). The total scat- 
tering rate for loss cone stars and the loss cone refilling rate have 
been estimated by various authors (e.g. Frank & Rees 1976; Light- 
man & Shapiro 1977; Rees 1988; Magorrian & Tremaine 1999; 
Merritt & Poon 2004; Wang & Merritt 2004; Ivanov, Polnarev & 



Saha 2005; Holley-Bockelmann & Sigurdsson 2006) and may vary 
over several orders of magnitude depending on the assumptions 
used. We use the scattering rate of Wang & Merritt (2004), 

which utilises the M-ct relation of Merritt & Ferrarese (2001), and 
assume that M BH is 3.61 x 10 6 M S (Eisenhauer et al. 2005). Whilst 
the above rate is derived for solar-type stars only, it scales with mass 
and radius only as m, l ^r\^ . For the models we have run we find 
that the overall scattering rate integrated over realistic populations 
differs from this by less than a factor of 3, provided the slope of the 
upper IMF is Salpeter-like. The amount of difference depends on 
the IMF, with more top-heavy IMFs having smaller scattering rates. 
In the case of IMFs with flatter upper slopes, such as that posited 
for the core of the Arches cluster (Stolte et al. 2005), the difference 
may be a factor of 10 or more. 

We assume stars, initially distributed according to either a 
Salpeter or top-heavy IMF, are supplied to the loss cone at a point 
in their life chosen randomly but weighted according to the mass 
and radius distribution above. This results in a number of distinc- 
tive characteristics for those stars which end up on loss cone orbits. 
Firstly, short lifetimes and the mass-dependence of the scattering 
rate ensure that few of the most massive stars (> 15M Q ) contribute 
to the population, even for top-heavy IMFs. Since such stars, once 
stripped, would probably display WR-like emission lines in their 
spectra unlike the observed Sgr A* cluster stars, this seems reason- 
able. Secondly, if one assumes a Salpeter IMF, most of the loss cone 
stars are in fact low mass main sequence stars which will be dis- 
rupted completely by a tidal-radius encounter with the black hole. 
Even if one assumes an IMF weighted towards more massive stars, 
as might be typical of the Galactic Centre region, most of the stars 
which are stripped will be main-sequence stars on close orbits. The 
Sgr A* cluster stars in this picture therefore arise from a relatively 
small proportion (a few percent) of the total number of stars scat- 
tered; those on the giant or horizontal branches. 

A further complication arises in that tidal stripping is not the 
only possible outcome of stars passing close to the black hole. Stars 
which pass the black hole at a distance above a tidal radius may 
still dissipate energy due to tidal distortion of the envelope, becom- 
ing tidally captured but not stripped. Such interactions may occur 
for initial peribothron distances of about three tidal radii (Novikov, 
Pethick & Polnarev 1992). On subsequent peribothra further tidal 
interactions take place, gradually causing the orbit to decay. If the 
star were able to survive these repeated interactions (which is un- 
likely, see Novikov et al. 1992) this would eventually result in a 
very close, circular orbit around the black hole. Relatedly, some 
main sequence stars are also scattered onto the red giant loss cone. 
On their first passage past the black hole these stars do not come 
close enough for significant tidal interaction at their current radii 
but they subsequently evolve to large enough radii that stripping 
can occur on a later orbit (Syer & Ulmer 1999). It is possible that at 
least some of the loss cone orbits are chaotic and indeed that chaotic 
orbits may be required if the loss cone is to be kept full (e.g. Merritt 
& Poon 2004, Holley-Bockelmann & Sigurdsson 2006), in which 
case this population may be minimal. However, if not, the popula- 
tion of stars late-stripped in this way should be larger than the pop- 
ulation from direct stripping on first peribothron passage. Since the 
radius-dependence of the scattering probability is relatively shallow 
the magnitude of this late-stripped stars effect is less than might at 
first be assumed. However, it is important in that nearly all the stars 
supplied to the black hole in this way are giants which are much 
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Figure 2. Energy required to strip 90 % (solid lines) and 99 % (dashed lines) of a star's envelope over the later phases of its lifetime, for different initial 
masses. Also shown greyed-out is the area in energy corresponding to the range of orbital energies which Sgr A* cluster stars with known orbits (Eisenhauer 
et al. 2005) would have if they were the mass of the resulting stripped stars. The darker grey area indicates the region of overlap between the energy range for 
the 90 % and 99 % stripped cores. 



less likely to suffer complete disruption than the main-sequence 
stars. Therefore stars which are stripped on their first pass of the 
black hole give very much a lower limit to the number of poten- 
tially observable Sgr A* cluster stars. An estimate of the popula- 
tion from this later stripping comes from running our simulations 
with the assumption that a star is stripped if scattered onto the loss 
cone corresponding to its maximum rather than current radius, i.e. 

^scatter 01 tn* '^'max,*- We then assign it a probability-weighted pe- 
riastron distance in the range r tidal!Cum , nt to r tidal!max . radius and add 
it into the stripped stars population when it is large enough to over- 
flow at periastron. This approach overestimates the population only 
in cases where the evolutionary timescale is shorter than the orbital 
timescale, which is not the case for most stars. It should be noted 
that we restrict the maximum-radius calculation to ages of less than 
10 10 years, so as to exclude stars which cannot have evolved to their 
present state in the time available, e.g. 0.5M e giants. This leads ef- 
fectively to a mildly top-heavy IMF with a turn-over at around 1M© 
amongst the stars which end up being stripped when a Salpeter-like 
IMF input population is assumed, since the maximum radii of more 
massive stars within the time limit are larger. 

Given the r-dependence for stars passing within their tidal ra- 
dius above and assuming stars typically increase their radii by a fac- 
tor of a few hundred between the main sequence and giant phases, 
it seems likely that the late-stripping giant population is roughly a 



few times the population of main-sequence stars stripped on their 
first pass of the black hole. This is borne out by simulation results; 
the factor by which the scattered population increases varies from 
less than 0.1 for stars with a 'normal' IMF (Kroupa 2001; the small 
change is due to the high proportion of low-mass stars) to over 3 
for stars with the most top-heavy IMF. However the top-heavy IMF 
populations had the lowest scattering rates to begin with. We call 
models incorporating stripping after the first pass 'late-stripping' 
models. As well as significantly increasing the visible Sgr A* clus- 
ter population and the number of non- visible remnants on similar 
orbits to the Sgr A* cluster stars in comparison to equivalent first- 
pass stripping models, adding in late-stripped giants can increase 
the accretion rate onto the black hole from tidal stripping by a factor 
of a few. Since there can be a significant time delay associated with 
lower-mass stars reaching their maximum radii, the late-stripping 
Sgr A* cluster population is likely to change with age and will ini- 
tially be much more top-heavy in mass than the input population. 



3.2 IMF 

Since the scattering probability varies inversely with mass, scatter- 
ing from an initially Salpeter-IMF reservoir of stars leads to a sig- 
nificantly bottom-heavy mass distribution for stars that pass within 
a tidal radius of the black hole. This is compounded by the short 
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lifetimes of the more massive stars. However, it is those more mas- 
sive stars which have sufficient luminosity to be observed in the 
GC region. Of a 'normal' Salpeter-like first-pass stripped popula- 
tion only around one star per 38000 is initially above 5M Q . With 
a scattering rate of ~ 5 x 10~ 4 stars yr -1 this obviously produces a 
very poor fit to the observed population of around 20 luminous Sgr 
A* cluster stars 1 . 

Though there remains the possibility that the scattering rate 
could be higher than the Wang & Merritt one (e.g. Holley- 
Bockelmann & Sigurdsson 2006), a large enough scattering rate 
would also violate other constraints imposed by observations of the 
GC region, requiring in particular an extremely large population of 
non- visible stellar remnants (Goodman & Paczynski 2005; see also 
section 4.4.1). 

There are two plausible effects which may help in increasing 
the Sgr A* cluster population without increasing the scattering rate, 
however. Firstly, many of the scattered stars come from the cen- 
tral region of the Galaxy itself. The simulations of Merritt & Szell 
(2006) suggest a lower limit of 10~ 4 stellar tidal disruptions per 
year for stars from the central cusp region, i.e. the inner few par- 
sees. There are numerous indications that the GC region may have a 
top-heavy IMF (Nayakshin & Sunyaev 2005). The present-day MF 
of the Arches cluster (around 30pc in projection from the central 
black hole) may have a turn-over at masses as large as 6M Q (Stolte 
et al. 2005). Top-heavy IMFs with lower mass limits of 1M or 
more can be a natural consequence of star formation in dense gas 
with a high opacity, such as in the GC region (Larson 2006). 

Secondly, as discussed above, adding in a population of late- 
scattered stars favours more massive, faster-evolving stars which 
reach larger radii. 

As well as running our simulations with a 'normal' IMF (we 
use the IMF given in Kroupa (2001), which has a Salpeter-like 
slope for massive stars and a lower limit of 0.01M Q ), we also run 
simulations with input IMFs which have lower mass limits at 1M Q 
(leading to one in 2500 stars above 5M Q ), 3M Q (one in 7 stars 
above 5M Q ) and 6M Q , and a hybrid IMF which assumes a rate of 
10~ 4 disruptions yr~' from a GC IMF with a turnover at 6M© and 
the remaining 4 x 10~ 4 disruptionsyr~' from a normal IMF. For 
the 3M e lower limit IMF we also run models with a relatively flat 
slope for stars above 6M , as is posited for the Arches cluster. 

Whether the top-heavy IMF models are realistic or not de- 
pends on the relaxation time in the central region, since these mas- 
sive stars still have to get on to the loss cone. The lives of more 
massive stars are short. If the relaxation time is much greater than 
the main-sequence lifetime of potential Sgr A* cluster-star progen- 
itors, then it is unlikely that enough of them are scattered. In partic- 
ular, if the relaxation time is 10 9 years then this suggests very little 
tidal stripping is possible from our top-heavy IMFs. Realistically, 
we require either a relaxation time below 10 8 years or else that the 
current spatial distribution of stars in the central region (in particu- 
lar the two inclined disks rotating in opposite directions) has some 
effect on the rate of stars scattering towards the black hole. The high 
scattering rate from the cusp region calculated by Merritt & Szell 



1 A 5 - 6M.0 star stripped on the RGB can live for 1 - 2 x 10 7 years; since 
the scattering probability increases with radius, about a third of the scattered 
> 5M stars are post-main sequence. For about 5 x 10 s - 10 6 years of this 
time they are observable in the K-band under conditions appropriate for the 
Galactic Centre. As stars below this mass spend shorter periods of their lives 
visible and blue, a rough calculation gives an Sgr A* cluster population over 
300 times too small. 
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Figure 3. HR diagram showing K-band contours (dashed lines) from K = 10 
to 20. The observed Sgr A* cluster stars lie between 14 and 16.8 (Eisenhauer 
et al. 2005). Also shown are unstripped stellar tracks for 2.7, 4, 5.9, 8.8 and 
13.3 M stars (grey) and, for reference, a portion of the stripped track for 
the 8.8 M RGB base star shown in figure 1 (black solid line). 

only includes the contribution of stars at r < 0.7 pc, so the relax- 
ation time of this region is appropriate here. Schodel et al. (2003) 
suggest the high stellar density in the central arcsecond should lead 
to a relaxation time of less than 10 8 years, and this value may be 
reduced further if one takes into account the potentially-large pop- 
ulation of compact remnants which are likely to also exist there 2 
(Davies & King 2005) and the possible effects of resonant relax- 
ation (Rauch & Tremaine 1996; Hopman & Alexander 2006). A 
lower limit slightly further out is provided by the unrelaxed stellar 
disks at ~ 0.4 pc. However, the presence of Wolf-Rayet stars and 
LBVs suggests they are only a few 10 6 years old. The top-heavy 
IMF scenarios are therefore at least not ruled out. If on the other 
hand they are not possible, then the alternative (given the evidence 
for a top-heavy IMF in the central region) is not lower-mass star 
scattering from the central region but a reduced overall rate of scat- 
tering in which all the loss cone stars come from much further-away 
regions with more normal IMFs. 



4 OBSERVABILITY 

Observations of stars in the Galactic Centre are only possible in 
the infrared due to the heavy visual extinction (some 35 magni- 
tudes) of the area. This presents particular problems for hot, blue 

2 Note though that if the density of objects in this region is too high then re- 
peated interactions could strip some of the envelope from input stars before 
their interaction with the black hole. This would reduce the orbital energy 
budget and skew the resulting orbit distribution towards higher eccentrici- 
ties. 
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stars for which most of the energy output is far from the observ- 
able region of the spectrum. Nevertheless, H and K-band spectra 
have been obtained for some of the Sgr A* cluster stars (Ghez et al. 
2003, Eisenhauer et al. 2005). For the brightest, SO-2 (S2 in Eisen- 
hauer at el. 2005), Ghez et al. (2003) observe the H I (4 - 7) Bry 
line at 2.166/jm and the He I triplet at 2.11^m; other lines in the 
region, such as N III at 2.1 155,um, are not seen. They obtain an ef- 
fective temperature for the star of around 30, 000 K. Using the hot 
star spectral atlas of Hanson, Conti & Rieke (1996), this led them 
to conclude that SO-2 was a late-type O star. Further observations 
by Eisenhauer et al. (2005) confirmed that the majority of the Sgr 
A* cluster stars have B star-like spectra. Interestingly, they appear 
to fall into strong-lined and weak-lined groups, suggesting some 
form of dual origin or dual evolution is going on. 

However, for many types of massive star there is a significant 
amount of overlap in IR spectral morphology (Morris et al. 1996) 
leading to potential confusions in classification. Whilst it is difficult 
to find template spectra for stars which have undergone similar evo- 
lution to our putative Sgr A* cluster-forming process, to a first ap- 
proximation they will appear similar to B stars (R. Napiwiotzki, pri- 
vate communication). Those stars which fall into the observable re- 
gion shown above have surface gravities, luminosities and effective 
temperatures in ranges which are similar to B stars. One important 
constraint imposed by the spectra is that a significant surface hy- 
drogen abundance must remain. However, whilst the stripped stars 
which remain luminous enough to be observed display He-enriched 
surface abundances after envelope stripping, nearly all retain a sur- 
face H mass fraction of between 0.1 and 0.6. Those which retain 
more envelope have higher surface H abundances; stars which lose 
their entire envelopes have, of course, rather low surface H (< 0. 12 
by mass fraction). 

To gain a rough idea of which stars will be visible, we cal- 
culate K-band limits assuming a distance modulus to the Galactic 
Centre of 14.6. For the K-band extinction we use the average value 
found by Blum, Sellgren & DePoy (1996), A K ~ 3.3. However the 
authors note that A K is rather variable, and can reach ~ 6 in places. 
Tabulated bolometric corrections and V — K colours are taken from 
Johnson (1966). The resulting K contours are shown in Fig. 3. It is 
apparent that red stars in the galactic centre region are significantly 
more observable than blue stars. 

In Fig. 4 we show the effect of IMF and amount of envelope 
remaining on the HR diagram positions of the stripped stars. It is 
apparent, as noted previously, that the 'normal' IMF populations 
cannot reproduce a luminous-enough group of stars. In addition, 
unless nearly all of the envelope is removed, they produce a pop- 
ulation of relatively low-luminosity red stars which will be visible 
in the K-band. It is thus unlikely that a normal-IMF first-pass pop- 
ulation will be able to reproduce the Sgr A* cluster stars even if 
the scattering rate is higher than previously assumed. However, the 
populations with top-heavy IMFs produce a number of luminous 
blue stars, although many of these populations also produce some 
red stars as well. For simulations including only first-pass strip- 
ping the visible populations are generally blue, but insufficiently 
numerous. This can be seen in snapshots taken at points during the 
simulations, as in Fig. 5. 

In Fig. 6 we show some output models from the late-stripping 
population, i.e. the population including main-sequence stars which 
scatter onto the giant loss cone and are stripped on a later pass of the 
black hole when they have evolved to larger radii. It is notable that 
the population is much-increased; the effective scattering rate for 
giant stars is much higher in this scenario. Models which are top- 
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Figure 4. Density distributions of overall Sgr A* cluster star HR diagram 
positions for varying assumptions about envelope stripping and IMF. Pan- 
els a and b show populations with 1% of envelope remaining and top-heavy 
(turnover at 6M Q ) and 'normal' Salpeter-like IMFs respectively. Panels c 
and d show the Salpeter IMF with 5% and 0% remaining envelope. Shown 
also are the approximate equivalents of the upper and lower K-band lumi- 
nosity limits observed for the Sgr A* cluster stars. 

heavy and which have less than 5 % remaining envelope produce 
reasonable Sgr A* cluster populations. 

4.1 Red Stars 

Some of the models which produce an acceptable population of 
blue stars also produce a small number of red ones. Whilst exclud- 
ing models which produce a few red stars does not rule every model 
out, it does rather restrict the parameter space over which the mod- 
els work. Of course, we could simply seeing the Galactic Centre 
at an unusual time when the red star population is low. However, 
the existence and observability of these model red stars may not be 
guaranteed. 

One possible reason that we might not see any red stars in the 
Galactic Centre, even though we predict some on close orbits, is 
that most such stars are convective and have relatively large radii. 
Even if they are not overflowing their tidal lobes, the magnitude of 
tidal effects on stars with convective envelopes is much greater than 
that on stars with radiative envelopes (see e.g. Hurley et al. 2002). 
This arises from the different dissipation mechanisms operating in 
each case. Stars which are not overflowing their tidal lobes but pass 
within a few r^l of the black hole are still subject to smaller-scale 
tidal interaction at periastron (Novikov et al. 1992). The likely out- 
come is a decrease in the eccentricity of the orbit, with the orbital 
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Figure 5. Snapshots of the HR diagram positions of some of our stripped 
star populations. Shown are population snapshots at age 4 x 10'yr assuming 
stripping leaves 1% of the stellar envelope remaining, with input IMFs with 
lower mass limits of (a) O.OIMq , (b) 3M with flat IMF slope, (c) hybrid 
IMF with 4/5 'normal' IMF, 1/5 3M e lower limit with normal IMF slope 
and (d) 6M lower limit. Also shown are the approximate luminosity limits 
corresponding to the K-band magnitudes for observed Sgr A* cluster stars, 
as detailed above; our populations should ideally contain ~ 20 stars in the 
blue region of this area (S2 has T e ft ~ 30000X', Ghez et al. 2003), not too 
many above the upper line (though the variable extinction may allow the 
upper line to be increased by a few magnitudes), and no red stars above 
either line. The 'normal' IMFs fail to produce anything like a Sgr A* cluster 
star, whereas the other IMFs produce acceptable distributions (some blue 
stars, no red stars) but at too low a rate. 

energy thus liberated deposited in the internal energy of the star. 
This in turn leads to a radius increase and hence stronger interac- 
tion on the next periastron passage. Novikov et al. suggest that this 
could lead to the destruction of the star in a few orbits, with the 
energy deposited exceeding the total binding energy. Such a mech- 
anism could efficiently rid the Sgr A* cluster of its red members. 
For example, if we remove from our sample of stars those which 
have convective envelopes and pass within 3 r t ^\ at pericentre, this 
effectively gets rid of all but the least massive red stars. 

A further complication is the effect of their extreme envi- 
ronment on their appearance. It has been suggested (Jimenez et 
al 2006) that irradiation of giant stars by even a relatively low- 
luminosity AGN can, via the destruction of molecules in their at- 
mospheres, make them appear significantly bluer. They claim such 
stars would have a similar appearance to the weak-lined Sgr A* 
cluster population. Such a claim should be handled cautiously, 
since irradiation typically affects predominantly the side of the star 
facing the radiation source, and it is unlikely that all the weak-lined 
Sgr A* cluster stars are fortuitously presenting their irradiated faces 





Figure 6. As Fig. 5, but showing successful and unsuccessful late-stripping 
populations: (a) Top-heavy (6M e limit) IMF with 1 % remaining envelope; 
(b) the same but with % remaining envelope; (c) Top-heavy (3M limit) 
IMF with 5 % remaining envelope, which produces a too-red population; 
(d) 'normal' IMF with 1 % remaining envelope, which again produces a 
large red population. 



in our direction at one moment. It does, however, suggest that even 
if we produce what should look like giants on Sgr A* cluster star 
orbits they are not guaranteed to appear to be so. Putting into our 
simulations the level of apparent surface temperature increase sug- 
gested by the simulations of Jimenez et al. we find a relatively mild 
effect on the overall populations, certainly not enough to pass off 
most red stars as members of the Sgr A* cluster. 

Overall, we find a significantly better fit with top-heavy IMF 
populations and with smaller amounts of remaining envelope. In 
particular, if more than a small proportion of input stars are of rel- 
atively low mass we find too great a population of low-luminosity 
but still visible red stars, which cannot be explained away by either 
of the above mechanisms. In fact, the absence of such observed 
stars could place an upper limit on the scattering rate of stars from 
outside the central region (see section 5). Whilst the orbits that we 
find for these low-mass stars are wider and more eccentric than both 
those of the more massive stars and the observed Sgr A* cluster or- 
bits, implying that they spend most of their lives far from pericen- 
tre where they may be less observable, we generate a large enough 
number of them in these cases that we would still expect some to 
be visible. That we do not suggests that, if this model for the Sgr 
A* cluster is to work, the input population is largely drawn from the 
central region itself. For example, late- scattering with a 6M e lower 
IMF limit produces a plausible Sgr A* cluster star population (Fig. 
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7) for the scattering rate we use, and in general input IMFs which 
are lacking in stars below 3M Q can be made to produce plausible 
Sgr A* cluster populations for scattering rates which are not much 
different. 

The division by Eisenhauer et al. (2005) and others of the Sgr 
A* cluster star spectra into two classes suggests the possibility of 
two separate observable populations. There are a number of possi- 
ble explanations for this within the framework of our model. Firstly, 
it is interesting to note that the simulations of Jimenez et al. (2006) 
suggest spectra similar to those of the weak-lined population may 
be produced by the irradiation of giants. In some of our simulations 
we find some giants amongst our Sgr A* cluster stars. It could be 
that, when irradiated, these make up the weak-lined population and 
the initially-blue stars make up the strong-lined population. How- 
ever as noted previously, the sort of temperature increases they pre- 
dict are not enough to make a significant difference to the red stars 
that we find. It is also possible, if some of the input stars are sup- 
plied by a top-heavy IMF and some are supplied by a 'normal' 
IMF, that we are seeing two populations with different input mass 
distributions. Finally, the evolutionary tracks of many of our more 
luminous stripped stars take them into the K-band visible region on 
the HR diagram twice: once moving bluewards shortly after strip- 
ping, and once moving back a little way redwards after a period on 
the helium main sequence. 

4.2 Orbits 

For both potential populations, the resulting post-stripping orbits 
are governed by the fraction of the envelope which is lost in tidal 
disruption, as this determines the energy budget available. Detailed 
SPH simulations are ongoing to determine the outcome of the strip- 
ping process (Dale et al., in preparation). At present it is an un- 
known, so we treat the amount of envelope remaining as a free 
parameter. Constraints are introduced, as in the previous section, 
by the requirement that the stripped star population is luminous 
enough to be observable (difficult if no envelope is left) and mostly 
blue (difficult if 5% or more of the envelope is left). This may be 
seen in Fig. 4. 

We consider primarily the situations in which 5, 1 and 0% of 
the envelope remains. Since a great deal of the envelope binding 
energy is concentrated in the region nearest the core, these sce- 
narios have significant differences in available energy budget. In 
particular, for stars with 5 % or more of the envelope remaining, 
the relatively low amount of energy extracted leads to a distribu- 
tion of orbits which are rather weakly-bound. As the periastron dis- 
tance is constrained by the need to initially pass close to the central 
black hole, this leads to very high typical eccentricities - 0.99 or 
more. The same applies to populations with a large proportion of 
low-mass stars, which also have smaller envelope binding energies. 
This includes the 'normal'-IMF population. These highly-eccentric 
orbits provide a poor match to the sample of Sgr A* cluster orbits 
which are known (see e.g. table 2 of Eisenhauer et al. 2005). Whilst 
the Sgr A* cluster orbits are eccentric, their measured eccentrici- 
ties range from 0.36 to 0.94. However, the models which produce 
unrealistic orbits are also those which have been previously ruled 
out for producing too many red stars. Stars from top-heavy IMFs 
and those which have more envelope removed, favoured as Sgr A* 
cluster star progenitors in the previous section, also provide the best 
match here. The greater the amount of envelope removed, the better 
the fit we find to observed orbits. Shown in Fig. 8 are snapshot or- 
bital distributions for some of our models which produce plausible 
luminosity distributions. 



Although the parameters of three of the six Sgr A* cluster stars 
which have detailed orbital solutions (S12, S14 and S8) are easy to 
match with observable models, the star with the best-determined or- 
bit, S2, is slightly outside the area in eccentricity and periastron oc- 
cupied by our models, and the two relatively low-eccentricity stars, 
SI and S13, are well outside this region. This suggests later orbital 
evolution after the initial stripping event. As discussed in Davies & 
King (2005), such orbits could be produced by scattering interac- 
tions with other stars in the central region after the initial stripping 
event - a not entirely surprising happening since, as we have al- 
ready discussed, the relaxation time in the central region needs to 
be small for the model to work in the first place. The relaxation time 
in the central arcsecond is likely to be less than 10 8 yr (Schodel et 
al. 2003), perhaps as low as 10 7 yr (Davies & King 2005) given the 
large number of compact remnants (white dwarfs, plus some neu- 
tron stars and black holes) which are likely to be there. Typically, 
many of the cluster stars which are observable in the top-heavy IMF 
model live around 1 - 2 x 10 yr after stripping. The portion of the 
tracks for which the star is visible varies from star to star but, whilst 
some stars are visible for a period of time just after stripping, a sub- 
stantial portion of the visible time generally occurs near the end of 
the stripped star's life. Thus it is plausible that some of the Sgr A* 
cluster stars have undergone one post-stripping scattering event. 

Another possibility is that, even though the Sgr A* cluster stars 
do not overflow their tidal lobes at pericentre, smaller-scale tidal 
effects are still occurring, leading to a slow circularisation of the 
orbit similar to orbital circularisation in massive binaries. 

4.3 Rotation 

The Sgr A* cluster stars have been observed to have rotation rates 
typical of main-sequence B stars (Eisenhauer et al. 2005). In our 
scenario for their origin, the rotation rate depends largely on the 
details of the envelope removal process. Since we are here treating 
that process as a 'black box', it is not possible to comment too 
deeply on the rotational rates. However, it is likely that the stripping 
process will induce at least some rotational velocity change in the 
star. Diener et al. (1995) find the effect of tidal interaction on their 
initially non-rotating model stars is to produce both rotation and 
oscillation. 

4.4 Further Population constraints 

Besides producing Sgr A* cluster stars like those observed, there 
are a number of other common-sense constraints on which popula- 
tions are acceptable. 

4.4.1 There should not be too massive a collection of remnants 
orbiting the black hole 

One of the chief complaints of Goodman & Paczynski (2005) re- 
garding an AGB star origin for the Sgr A* cluster stars is that a 
sufficiently large population of stripped AGB stars (which only re- 
main in the observable region for a short time) to reproduce ob- 
servations would lead to an unreasonably large mass of non-visible 
remnants orbiting the black hole (white dwarfs, neutron stars and 
some black holes, depending on the IMF). Whilst the stripped RGB 
and horizontal branch stars which make up the bulk of the Sgr A* 
cluster in our scenario remain visible for a longer time, this is still 
potentially a problem. Several factors affect the mass of remnants 
in the central region. Firstly, the amount of mass which enters the 
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region has a strong effect. Top-heavy IMF models generally leave 
more massive remnants, so even if the number of stars coming in 
is the same the total remnant mass is greater. For example, 'nor- 
mal' first-pass IMF models leave a few hundred to a few thousand 
solar masses of orbiting remnants, generally towards the lower end 
of this range. This is because most stars which encounter the black 
hole are disrupted entirely. For top-heavy IMFs the remnant mass 
can be closer to a few 10 5 solar masses. However the total number 
of remnants is similar to that for 'normal' IMFs; the difference is 
that each remnant is more massive. 

Second, the mass of orbiting remnants is greater for the late- 
stripping populations. This is unsurprising, since a much larger pro- 
portion of them encounter the black hole as giants and so avoid 
complete shredding. For late-stripping top-heavy IMF populations 
the remnant mass resulting can be particularly unrealistic and can 
approach the mass of Sgr A* itself in the most extreme case. This 
is problematic, since the late-stripping top-heavy IMF models pro- 
vide the best match to the K luminosity and orbital characteristics. 
However, the remnant mass is dependent on the scattering rate. For 
Salpeter- slope IMFs with a high lower mass limit, the scattering 
rate is not very different from the 1Mq,1Rq scattering rate; for 
IMFs which are top-heavy because they have a flatter upper slope, 
such as our Arches-like IMF, the scattering rate is much lower and 
acceptable populations can still be obtained with reasonably small 
remnant masses. In addition, the above analysis assumes the loss 
cone is kept full - probably an unlikely state. If instead the scat- 
tering rate is constrained by the refilling rate of the loss cone the 
rate of disruptions may be much lower and hence the remnant mass 
still reasonable, although this will also reduce the number of visible 



stars unless formation happens in bursts, perhaps related to recent 
star formation in the GC region. However, this latter scenario would 
have trouble explaining the current low luminosity of Sgr A*. 

More importantly, the above analysis of the remnant mass as- 
sumes nothing happens on the formation of a remnant and no fur- 
ther interaction between remnants takes place after formation. In 
the top-heavy IMF case, much of the mass of remnants is made 
up out of neutron stars and even some black holes. This requires 
type II (potentially lib or even lb, since only a small amount of 
envelope remains) supernovae. It is generally thought (e.g. Lyne & 
Lorimer 1994) that the formation of a neutron star or black hole by 
a supernova is accompanied by a large velocity kick arising from 
asymmetries either in the mass ejection or the neutrino emission. 
Therefore much of the remnant mass could in fact have been ejected 
from the central region, either inward (i.e. the remnant is propelled 
onto an unstable orbit, from whence it is swallowed) or outward 3 . 
If all such remnants are ejected, the total mass of orbiting remnants 
in the top-heavy late-stripping models after 10 years decreases to 
~ 10 5 M Q . For non-expelled remnants, interactions arising from the 
high density of the region may also lead to ejections. In particular, 
feedback interactions may increase the number of ejections, with 
temporarily-overdense regions having high ejection rates; and, if 
there is a population of massive black holes which have migrated to 
the central parsec by dynamical friction (Miralda-Escude & Gould 



Note that the velocity dispersion in the central region is rather large; or- 
bital speeds may be > 1 000km s _1 . The average kick value from the study 
of Lyne & Lorimer is 450km s" 1 . So this mechanism probably accounts for 
only a portion of the remnants. 
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Figure 8. Snapshot of Sgr A* cluster orbits at 4 x 10 9 years, for (a) late-stripping IMF with Salpeter slope and lower limit of 6M , 1 % remaining envelope; 
(b) the same but with % remaining envelope; (c) Arches-type IMF with 1 % remaining envelope; (d) First-pass stripping with lower IMF limit of 3M and 
1% remaining envelope. Black points indicate visible, grey points non- visible Sgr A* cluster stars. Note that in general the visible stars are spread relatively 
evenly amongst the distribution of orbits. Shown also are the six current observed Sgr A* cluster orbits (Eisenhauer et al. 2005). 



2000) these will tend to eject the individually less-massive rem- 
nants from stripped-star evolution. If mergers between the central 
black hole and other Mailing SMBHs from Galaxy mergers hap- 
pen periodically (Yu & Tremaine 2003), the central region is likely 
to be 'cleared-out' of its population of remnants when this occurs: 
in this case, the central remnant population may be only the much 
smaller number of compact objects formed since the most recent 
SMBH merger. 

The number of remnants is also partially dependent on the 
amount of envelope which the stripping process typically removes 
from a star. The greater the amount of envelope the star retains, the 
more likely it is to undergo a second expansion phase after cen- 
tral helium burning and end up as a helium giant. In this state it is 
likely to overflow its tidal lobe a second time or at least undergo 
strong tidal interactions (section 3.1). Several outcomes are possi- 
ble for this second stripping process. If it behaves similarly to the 
first, one would expect a more circular orbit with the same pericen- 
tre distance to result. More probably, a second large input of energy 
and/or repeated tidal interactions at periastron may disrupt the star 
completely (Novikov et al. 1992). In either case, the total mass of 
remnants is reduced. 



4.4.2 The flaring rate should be reasonable 

The disruption of a star and subsequent accretion of its matter can 
produce a strong X-ray flare, such as have been observed in some 
AGN (e.g. Piro et al. 1988; Brandt, Pounds & Fink 1995; Grupe, 
Thomas & Leighly 1999; Donley et al. 2002). However, the time- 



scale of the tidal disruption event determines its observability, since 
very long timescale events are not detectable as flares. Red giants 
with large radii are amongst the stars which have a long disrup- 
tion timescale (Syer & Ulmer 1999). Therefore a larger disruption 
rate in which the increase is primarily of giant stars only, such as 
we find in the late-stripping case, does not necessarily result in a 
greater number of observable flares. Much of the X-ray luminosity 
function of active galactic nuclei may be powered by this method, 
particularly at the low-luminosity end (Milosavljevic, Merritt & Ho 
2006). Tidal disruption rates have frequently been calculated in the 
past with a view to working out the rates of X-ray flares (e.g. Rees 
1988; Magorrian & Tremaine 1999). More frequent minor flares 
have also been observed in Sgr A*, which could be the signatures 
of stars passing through the black hole's accretion disk (Nayakshin, 
Cuadra & Sunyaev 2004). 

Even for our most extreme late- scattering models, the rate 
of tidal disruption is never more than a factor of 4.5 greater than 
for the case in which all stars are 1M , 1R and are stripped on 
their first pass. This corresponds to a maximum rate for the Galac- 
tic Centre of around 2 x 10 _3 disruptionsyr _1 . For realistic mod- 
els which provide a reasonable fit to the observational data and 
to the other constraints discussed here the rates of tidal disrup- 
tion are generally lower. For example, the total rate of tidal dis- 
ruptions in the case of the 'Arches' IMF is only 20% of that for 
the IMq, 1R©, first-pass stripping case. However, as yet the con- 
straints on flaring rates are relatively weak. Donley et al. (2002) 
find the rates to be approximately 9 x 10~ 6 galaxy~' yr~' for inac- 
tive and 8.5 x 10~ 4 galaxy~ 1 yr~ 1 for active galaxies. Our predicted 
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flaring rates are nearly all comfortably in this range. This is there- 
fore also a rather weak constraint on the models when compared to 
the others. 

A further observational limit which has the potential to con- 
strain the rate of inspirals into the black hole in the future is their 
gravitational wave emission (e.g. Hopman & Alexander 2006). 

4.4.3 The accretion rate should be reasonable 

The tidal disruption process is also a method of black hole feed- 
ing, and it is likely that at least some fraction of the stripped matter 
(perhaps around half; Rees 1988) is accreted onto the black hole. 
Although the current low luminosity of Sgr A* indicates that its 
accretion rate at the moment is much smaller than expected (e.g. 
Baganoff et al. 2003; Zhao, Bower & Goss 2001), it must obvi- 
ously have accreted enough matter in the past in some fashion to 
grow to its current mass. Tidal shredding is by no means the only 
method of feeding a black hole; at least some of the mass is likely 
to arise from gas accretion (e.g. Wada 2004). If the amount of mat- 
ter accreted over 10 10 years from tidal disruptions is significantly 
over 3.6 x 10 6 Me, this suggests the model is unrealistic. Of course, 
the amount of accretable matter from loss cone stars changes over 
the lifetime of the black hole, as the scattering rate depends on 
the black hole mass. Wang & Merritt (2004) find, perhaps coun- 
terintuitively, a wider loss cone and thus greater scattering rate 
onto lower-mass black holes. The fraction of matter which is ac- 
creted is also a matter of debate and may be rather lower than half, 
perhaps closer to 10% (Ayal, Livio & Piran 2000). Late-stripping 
models in general have larger amounts of accreted mass simply 
because more stars are stripped. Similarly top-heavy IMF models 
have larger amounts of accreted mass because the total amount of 
mass in the stars encountering the black hole is greater. Generally 
the models which can be rejected on these grounds are the same 
ones which can be rejected on the grounds of having too great a 
mass of remnants, therefore. Top-heavy late-stripping models can 
produce produce acceptable levels of accretion if the scattering rate 
is not too high (this occurs with flatter upper IMF slopes, for exam- 
ple, so the Arches-like model is again good here) or if the fraction 
of matter accreted is closer to the Ayal et al. (2000) value than to a 
half. 



5 DISCUSSION 

If the Sgr A* cluster stars arise from relatively low mass progeni- 
tors, the 'paradox of youth' (Ghez et al. 2003) is no longer a prob- 
lem. However, it seems apparent from our models that the require- 
ment for young stars is only slightly relieved; the model which 
matches most closely to the HR diagram position of a 15M Q ZAMS 
B star, for example, is an initially 9M.g, star which undergoes strip- 
ping at an age of 29.5 Myr to remove all but about 2M© of its 
mass. But the 15M e star itself stays on the main sequence for 10 
Myr. If the Sgr A* cluster stars are produced in this manner, then, 
it is likely that the disks of apparently massive young stars within a 
parsec of Sgr A* which are the next furthest objects out from them 
were in fact formed in place (Nayakshin & Sunyaev 2005) and that 
the nearby Galactic Centre region in general supplies many of the 
scattered stars which become visible on Sgr A* cluster orbits after 
stripping. 

The IMF of the scattered stars is the most important ingredi- 
ent in producing a plausible Sgr A* cluster model from stripped 
stars. Most top-heavy IMF models can be made to produce enough 
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Figure 9. Cumulative number of visible red stars with radius expected in the 
Galactic centre from 'normal'-IMF late-scattered stars after 4 x 10 9 years, 
with 1 % remaining envelope and tidal disruption rates (upper line to lower 
line) 1(T 3 , 5 x lCr 4 ,lCr 4 and 10" 4 disraptionsyr"' . 

visible blue stars and few enough red stars with scattering rates 
which are in the range of those claimed in the literature. However, 
as discussed above, other constraints are important. In particular, 
the fraction of the matter produced by stripping stars which is ac- 
creted onto the black hole should be relatively low - a quarter or 
less - if this model is not to produce a greater mass for Sgr A* than 
is observed, and the relaxation time in the central 0.7 pc needs to be 
below 10 8 yr. The most successful model when all these constraints 
are taken into account is the late-stripping Arches-like' model. If 
the constraints are not met, it is probable that this mechanism still 
operates at least at some level, and could therefore be responsible 
for some of the Sgr A* cluster stars. 

The generally poor results we find with 'normal' IMFs may 
also provide an interesting constraint on the tidal disruption rate 
from regions outside the central cusp. Stars scattered from these 
regions are likely to belong to such an IMF - in fact, it is likely 
that they are entirely deficient in high-mass stars due to the short 
lifetimes involved. Even though lower-mass stars end up on higher- 
eccentricity orbits, in the late-stripping case enough red stars are 
made that we would still expect a population of red stars observ- 
able amongst the Sgr A* cluster stars. The fact that we do not 
see such a population potentially suggests limitations either on the 
tidal stripping rate from a full loss cone or the rate stars are scat- 
tered onto the loss cone. In Fig. 8 we plot the number of visible 
red stars likely with radius as a function of tidal disruption rate. A 
significant red population arises from scattering rates above about 
10 _4 disruptionsyr _1 . 



5.1 Observational Prospects 

Observationally, there are a number of ways that it may be pos- 
sible in the future to distinguish between this scenario and others 
for the formation of the Sgr A* cluster. An unavoidable part of the 
stripped-star scenario is that there is a large population of single 
compact objects, many of which will be neutron stars or even black 
holes, in the central region. This extra dark mass should cause the 
orbits of the cluster stars to precess slightly (Mouawad et al. 2005; 
Weinberg, Milosavljevic & Ghez 2005). Whilst current telescopes 
are unable to measure the stellar orbits with enough accuracy to 
make out this effect, future 30 - 100 m Extremely Large Tele- 
scopes (ELTs) should be able to effectively weigh the extended 
matter distribution around the black hole using this method (Wein- 
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berg et al. 2005). There is also the possibility of observing these 
remnants directly. Pulsars orbiting Sgr A* should be visible with 
SKA (Cordes et al. 2005). Recent work (Wu et al. 2006; Deegan & 
Nayakshin 2006) also suggests that single neutron stars and black 
holes may be observable in the Galactic Centre region if they can 
accrete enough gas from the local ISM to be X-ray bright. This 
effect has been suggested as an explanation for the overdensity of 
X-ray point sources towards the Galactic centre (Muno et al. 2003). 
However, although the absence of a cusp of compact objects would 
disprove the stripped-stars theory, the presence of a such a cusp is 
not in itself definitive proof as massive compact objects (5 - 10M e 
black holes) are expected to migrate to the central parsec from out- 
side via dynamical friction (Miralda-Escude & Gould 2000), and 
such objects may be efficient at ejecting lower-mass stellar rem- 
nants from the area. 

With future space-based or AO telescopes performing deeper 
observations in the central parsec region, there exists the possibility 
of constraining the initial or at least current mass function in this 
area. As is apparent above, the stripped-stars model only works if 
the local IMF is top-heavy. In particular, the best-fit IMF is Arches- 
like', with a relatively flat upper slope. There are already limits on 
the size of any potential low-mass star population from X-ray emis- 
sion (Nayakshin & Sunyaev 2005), and the extremely high WR and 
LBV population in the central parsec in comparison to other mas- 
sive clusters with a similar O and B star population (e.g. Tanner et al 
2006) suggests that the upper IMF slope may be flatter than normal. 
In our simulations we find a distribution of cluster star luminosities 
which extends down below the boundary of what is currently ob- 
servable, so we would expect future deeper observations of the Sgr 
A* cluster itself to reveal more blue stars on similar orbits, as indi- 
cated in figure 7. 

Finally, there are a number of properties of the cluster stars 
themselves which could be used to confirm or deny different theo- 
ries as to their origin. Whilst the stripped cores we have produced 
appear similar to main-sequence B stars in some observational re- 
spects, they should be by no means identical. In particular, their 
surface compositions will differ. If the Sgr A* cluster stars are B 
stars then their surface H and He abundances should be typical of 
the metallicity of the Galactic centre region - i.e. about solar (e.g. 
Najarro et al. 2004). If they are stripped stars, then they should be 
somewhat H-depleted, probably with a surface mass fraction of hy- 
drogen between 0.6 and 0.1, and these surface abundances should 
vary strongly from star to star. As the orbit is dependent on the 
amount of envelope removed, there may be correlations between 
the periastron and eccentricity and the luminosity of the Sgr A* 
cluster stars which could be investigated in the larger sample of 
stars and orbital solutions available from deeper observations. Here 
we would expect a slight correlation between high luminosity and 
lowish eccentricity - however, this may be wiped out entirely if 
further scattering events take place to form the low-eccentricity or- 
bits (section 4.2). A much more definite test would be if we were 
lucky enough to observe a supernova amongst these stars, since we 
expect some of our top-heavy IMF models to undergo type lib su- 
pernovae. B stars on the same orbits probably avoid supernovae; 
once they evolve to larger radii, it seems probable that they are 
completely disrupted by tidal interactions with the hole. 

5.2 Central Star Clusters in Other Galaxies 

It is interesting to note that the mass of the central black hole is crit- 
ical in this picture in deciding whether a central cluster like the one 
observed can be formed and what its characteristics are. There have 



already been observations suggestive of central clusters of young 
stars in other galaxies (e.g. Lauer et al. 1998), and it is likely that if 
the stripped-stars mechanism operates in our own Galactic centre 
it also operates in others. This raises the possibility of observing 
clusters analogous to the Sgr A* one around other SMBHs, and 
the question of what properties such clusters would have. Firstly, 
the scattering rate depends on ni(, BH (where m(, BH is the black 
hole mass in units of 10 6 Ma) roughly as m^ B ^ 5 (Merritt & Wang 
2004), i.e. a black hole of a few 10 8 M Q should, discarding any 
other effects, have a Sgr A*-like cluster only one-third as numerous 
as the Galactic one. Secondly, the Schwarzschild radius is greater 
for more massive black holes. For black hole masses comparable 
to that of Sgr A* the Schwarzschild radius is only a few solar radii 
and most solar-type stars passing within a tidal radius of the black 
hole are disrupted rather that directly swallowed. However the ratio 
of Schwarzschild radius to tidal radius is approximately 



2.1x10"*^ 
Rt n 



where m t and r t are the mass and radius of the scattered star in 
solar units (e.g. Beloborodov et al. 1992, Ivanov & Chernyakova 
2006). Black holes of a few 10 8 M e will therefore swallow solar- 
type stars directly rather than tidally disrupting them and accreting 
some of the gas. Many of the visible Sgr A* cluster stars in our 
first-pass models arise from horizontal branch (HB) stars, as this 
is the longest phase during which the star has a well-defined core. 
As HB stars are also the most compact of the potential progenitor 
stars, they will be the first section of the progenitor population to be 
swallowed rather than stripped as black hole mass increases. Taking 
typical values of mass and radius for HB stars which are visible 
as Sgr A* -like cluster stars when stripped suggests that the black 
hole mass above which the Sgr A*-like cluster population drops 
dramatically should be around 10 10 Mq - i.e. this effect will be 
very small in most cases. However, the fate of lower-mass stars and 
thereby potentially the fraction of gas from scattered stars which 
accretes should be significantly different at the upper and lower 
ends of the known SMBH mass spectrum. In particular, whilst the 
scattering rate may be lower for more massive black holes (Wang 
& Merritt 2004) there is a mild compensatory effect in terms of 
extra mass accretion. The observation of clusters obeying the mass 
relation above in other galaxies would act as significant evidence in 
favour of this method of forming the Sgr A* cluster. 

The effect of black hole mass on scattering rate is also impor- 
tant if tidal stripping is to be considered as a general AGN fuelling 
mechanism (Milosavljec, Merritt & Ho 2006). As noted above, the 
amount of mass supplied in this way is certainly sufficient for such 
a process. The decrease in loss cone angle as the black hole mass 
increases suggests that in this picture it is the early life of the mas- 
sive black hole which is characterised by more rapid accretion and 
that the growth of supermassive black holes gradually slows down 
with time. However, the evolution of massive black holes is un- 
likely to be this simple. If galaxy mergers are a frequent event then 
there are likely to be phases in which the black holes at the centres 
of some galaxies are binaries. Such a situation can cause a large 
increase in the tidal disruption rate (Ivanov, Polnarev & Saha 2005) 
and will also have a strong disruptive effect on any pre-existing Sgr 
A* -like cluster, with a high ejection rate from the central region (Yu 
& Tremaine 2003). 
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